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We show that the quantum Jarzynski equality generalizes to "PT-symmetric quantum mechanics with 
unbroken VT -symmetry. In the regime of broken VT -symmetry the Jarzynski equality does not hold as 
also the CPT-norm is not preserved during the dynamics. These findings are illustrated for an experi¬ 
mentally relevant system - two coupled optical waveguides. It turns out that for these systems the phase 
transition between the regimes of unbroken and broken "PT-symmetry is thermodynamically inhibited as 
the irreversible work diverges at the critical point. 
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Over the last two decades VT -symmetric quantum me¬ 
chanics has established itself as an important area of mod¬ 
ern research. In its original formulation [ 1 ] VT -symmetric 
quantum mechanics was a rather mathematical theory with 
only loose connection to physical reality [2]. However, re¬ 
cent experimental progress has provided a realization of 
PT-symmetric systems, which can be understood as sys¬ 
tems with balanced loss and gain [2-6]. In particular, Riiter 
et al. observed VT- symmetry in an optical system [7], 
which consists of two coupled waveguides with modulated 
refraction indexes. 

Around the same time of the first description of VT- 
symmetric quantum mechanics [1] Jarzynski achieved 
a major breakthrough in thermodynamics of small sys¬ 
tems [8]. The Jarzynski equality, (exp (—/3W)) = 

exp (—/3A F), allows to determine the free energy differ¬ 
ence for an isothermal process from an ensemble of non- 
equilibrium realizations of the process. Here, /3 is the in¬ 
verse temperature, W is the non-equilibrium work, and 
AT is the free energy difference. The angular brackets, 
(•), denote an average over all possible work values. There¬ 
fore, modern thermodynamics places special interest on the 
distribution of work values, V(W), which has been stud¬ 
ied extensively, for instance, in classical systems [9-12] as 
well as for quantum systems [13-22], 

It was shown that the classical Jarzynski equality can 
be generalized to isolated quantum systems [23, 24], for 
which the thermodynamic work is determined by a two- 
time energy measurement [25, 26]. In this approach one 
commonly considers the following procedure: a quantum 
system is prepared in a thermal Gibbs-state; then a projec¬ 
tive energy measurement is performed, before the system 
evolves under an externally controlled Hamiltonian, H t \ 
the procedure is concluded by a second, projective energy 
measurement [27]. 

The present work addresses the following question: 
Does the quantum Jarzynski equality together with the two- 
time energy measurements generalize to VT -symmetric 
quantum mechanics? We will find that this is, indeed, the 
case for systems with unbroken VT- symmetry, whereas 
for broken VT- symmetry the Jarzynski equality does not 


hold. These findings will then be carefully analyzed and 
illustrated for the model describing the optics experiment 
by Riiter et al. [7], 

We will show that the phase transition between unbroken 
and broken regime has a clear signature in the behavior of 
the irreversible work, (JT irr ) = (W) — AT. This is in 
analogous to systems exhibiting quantum phase transitions 
[28-30]. However, while quantum phase transitions are 
thermodynamically allowed, the transition from unbroken 
to broken VT -symmetry is thermodynamically inhibited. 

Fundamentals of VT-symmetric quantum mechanics. 
We start by briefly reviewing the main properties of VT - 
symmetric quantum mechanics. Consider a quantum sys¬ 
tem with non-Hermitian, but VT -symmetric Hamiltonian 
H, i.e., [PT, H] = 0. Here, V is the space reflec¬ 
tion (parity) operator, and T is the time-reflection operator 
[1,3,31,32], 

V xV =—x and V pV =—p 

_ ( 1 ) 
TxT = x, T pT = — p and TiT = — i 

where x and p are position and momentum operator, re¬ 
spectively. Since T also changes the sign of the imaginary 
unit i, canonical commutation relations such as \x. p] = ih 
are invariant under VT. It has been seen [2, 32] that VT - 
symmetric Hamiltonians generally exhibit two parametric 
regimes: a regime of unbroken VT -symmetry in which all 
eigenvalues of H are real, and a regime of broken VT - 
symmetry for which the eigenvalue spectrum has real and 
complex parts. 

The major difference between Hermitian and VT - 
symmetric quantum mechanics is the definition of the inner 
product [31, 32]. For Hermitian Hamiltonians we have, 

(V’llV’2) = Vl "02, (2) 

where, as usual, f denotes conjugate transpose. This in¬ 
ner product, however, yields indefinite norms for the non- 
Hermitian, but VT -symmetric case [1, 3, 31, 32]. This can 
easily be seen explicitly, for instance, for the system dis¬ 
cussed in the second part of our analysis. Therefore, the 
inner product in VT -symmetric quantum mechanics is de¬ 
fined in terms of the metric operator C as [31, 32] 

= (CVTifi) ■ ■ (3) 
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In the unbroken regime C can be determined from [3 1 , 32], 

[C, H]= 0 and C 2 = I. (4) 

It is worth emphasizing that the time evolution induced 
by a time-independent Hamiltonian with unbroken VT - 
symmetry is unitary, and hence the norm is preserved, 

{A\A) CV r = 1 for a111 [ 31 > 32 ]- 
The metric tensor C reminds us of the charge conjugation 
operator from held theory [32]. However, in the present 
case C is not necessarily associated with an observable, but 
rather defines a physically consistent theory. 

Quantum work for unbroken VT -symmetry. In the fol¬ 
lowing, we are interested in processes that are induced by 
a time-dependent control parameter A t , with H t = H (A t ). 
We start by considering a driving protocol for which H t = 
H{X t ) is in the regime of unbroken VT -symmetry at all 
times. Note that a time-dependent Hamiltonian H t also in¬ 
duces a time-dependent metric C t as can be seen from its 
definition (4). It has been recently shown [33] that then the 
generalized Schrodinger equation reads, 

ihd t \ip t ) = (H t + A t ) \ip) , (5) 

where A t is a time-dependent gauge held, that is necessary 
to preserve normalization under the time-dependent metric 
C t . It is given by, 

7 T) 

At = -- WA 1 d t w t , (6) 

where Wj is the transpose of the metric C t , i.e., W t = C[. 

Commonly, the thermodynamic work done during a pro¬ 
cess of length r is determined by a two-time energy mea¬ 
surement [25, 26]: at initial time t = 0 a projective energy 
measurement is performed; then the system is allowed to 
evolve under the generalized time-dependent Schrodinger 
equation (5), before a second projective energy measure¬ 
ment is performed at t = t. For a single realization of this 
protocol the work is given by 

^|0 m )->| <t> n ) = E n (K) ~ E m (X 0 ), (7) 

where \(p m ) is the initial eigenstate with eigenenergy 
E m (X 0 ) an d | cf> n ) with E n ( A T ) denotes the hnal state. 
The distribution of work values is given by averaging over 
an ensemble of realizations of the process, V(W) = 
(<5 (W — ^| 0 m >-)-| 0 n )))» which can be written as [20, 34] 

V(W) = ^5{w- W\4, m) ^ M ) p (| <p m ) -+ |<M) • 

m,n 

( 8 ) 

In the last equation the symbol ^ accounts for discrete and 
continuous parts of the eigenvalue spectrum. Without loss 
of generality we assume here that the spectrum is discrete, 
but see also Ref. [20]. 

In Eq. ( 8 ) p(\<p m ) —>• | <j> n }) denotes the probability to 
observe a specific transition \4> m ) —>■ | <p n ). This probabil¬ 
ity is given by [20, 34], 

p (| m) -0 |n)) = tr {n„ U r U m p 0 n m U ] T ) , (9) 


where p 0 is the initial density operator of the system and 
U T is the unitary time evolution operator, 

U T = T> exp dt. ( H t + At)j ■ (10) 

Here, 7> is the time-ordering operator, and 1 [„ denotes the 
projector into the space spanned by the v\h eigenstate. For 
the sake of simplicity we further assume that all spectra 
are non-degenerate, for which we simply have Ii u = (f) v ■ 
{CfPT <p v ) [35]. Hence, Eq. (9) can be written as, 

P(\(f>m) -f Ifn)) = ( CrVTfn) ' (U T (/> m ) 

• {CfPTfm) ■ (P0<f>m) ■ {CrVTUrfm) ■ <j) n . 

Now, let us assume that the system under study was ini¬ 
tially prepared in a Gibbs state, namely we have 

Po4>m = [exp (-/3 E m )/Z 0 \ (j) m , (12) 

where Z {] = tr {exp (—3H {] ) } is the partition function. 
Then, we compute the average exponentiated work, 

(exp (—/3W)) = f d WV(W) exp (~PW) 

„ J (13) 

= ^ exp {—f3E n + f3E m ) p (\(j> m ) -f |< p n }) . 

m,n 

Substituting Eqs. (11) and (12) into Eq. (13) and us¬ 
ing the C 'P T - n o mi a I i z a t i o n of the initial eigenstate 0 rn , 

(C 0 VT (p m ) ■ 4>m = 1, we have 

(exp (-PW)) = (1/Z Q ) exp (~/3E n ) 

m,n (14) 

X (CrVTtpn) ■ ( U T (p m ) ‘ (C t VT U T (p m ) ■ (p n . 

We further employ the CV'T -symmetric partition of the 
identity, I = ip v ■ (C,'P‘T ip v ) , which is invariant under 

unitary evolution [32], and the C'PT -normalization of the 
final eigenstate, (p n . Hence, we obtain 

(exp (-fiW)) = Z T /Z 0 = exp (—/3A F) , (15) 

where F = —1/P In (Z) is the free energy. In conclusion, 
we have shown that the quantum Jarzynski equality (14) 
remains valid for quantum systems with unbroken VT - 
symmetry, and for which the time-evolution is described 
by the generalized Schrodinger equation (5). 

Regime of broken VT -symmetry. What remains is to 
check, whether Eq. (15) is also valid in the regime of bro¬ 
ken VT -symmetry. Analogous to the unbroken regime the 
initial energy eigenstate can always be chosen to be CVT- 
normalized. Therefore, the treatment of the broken regime 
is identical to the above discussion, if we replace the defi¬ 
nition of C in Eq. (4) by [36] 

{C,H} = 0 and C 2 = I. (16) 

Comparing Eq. (4) with Eq. (16) we observe that in the 
regime of unbroken VT- symmetry the metric operator C, 
commutes with H, whereas in the broken regime the met¬ 
ric operator anticommutes with the Hamiltonian H [37]. 
In addition, the time-evolution in the broken regime ceases 


tonian as, 
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FIG. 1. (color online) Irreversible work (JTi rr ) (27) for the linear 
protocol (26) and ao = 1 /2, that is a system starting in the regime 
of unbroken PT-symmetry. Parameters are/3 = 1, h= 1, r = 1, 
and k = 1. Inset illustrates the validity of the PT-symmetric 
Jarzynski equality (15). 


to be unitary as the energy spectrum contains a complex 
part. Therefore, even when including a gauge field in the 
generalized Schrodinger equation (5) to account for time- 
dependent metrics ( 6 ), norms are not preserved. In particu¬ 
lar, we have i’m-{C T VTi/j m ) / I, with t/’m = U T <j) m , 

if the time-evolution operator, U T , is not at least unital 
[34, 38-40]. A unital map is a trace preserving, completely 
positive map, under which the identity is preserved, and 
that can be written as a superposition of unitary maps [41]. 

In conclusion, it becomes apparent by inspecting 
Eq. (14) that the quantum Jarzynski equality together with 
the two-time energy measurement approach does not hold 
under broken VT -symmetry. 

VT-symmetric Jarzynski equality in optics. The re¬ 
mainder of this discussion is dedicated to a careful analy¬ 
sis and illustration of the above findings for an experimen¬ 
tally relevant example. In a recent experiment Riiter et al. 
showed that VT -symmetric quantum mechanics can be re¬ 
alized in optical set-ups [7] - two coupled waveguides, of 
which only one is optically pumped. The optical-field dy¬ 
namics is described by [7] 


i d z E i — — Ei — n E 2 , 

id z E 2 = — — E 2 — n Ei , 


(17) 


where E\ :2 are the field amplitudes in the waveguides, k 
is the coupling constant, and 7 is the gain coefficient due 
to the optical pumping. Identifying the spatial coordinate 
with a time variable, 2 = t/h, the dynamics of a VT - 
symmetric quantum system is “frozen” in the profile of the 
field amplitudes. In particular, we can identify the Hamil- 


H(a) = k 


ia — 1 \ 

—1 —iaj ’ 


(18) 


where we introduced the new parameter a = j/2k. Note 
that this Hamiltonian belongs to the class of two-level sys¬ 
tems discussed extensively in the literature [3, 31, 32]. 
With the parity operator, 


V = 



(19) 


and noting that T performs here only complex conjuga¬ 
tion one easily convinces oneself that H (a) (18) is, indeed, 
VT -symmetric. The eigenenergies are given by, 

e lj2 = ±kV 1 - a 2 , ( 20 ) 


from which we conclude the regime of unbroken VT- 
symmetry, a < 1 , and the broken regime is a > 1 . 

It is then a simple exercise to determine the eigenstates 
and the metric operator C. We obtain in the unbroken 
regime for the eigenstates. 


IA") = 


1 <t>T) = 


V 7 2 Vl - a 2 V arcsin ( a ) 


g— ^ arcsin(a) 


1 


V2V1 - a 2 

with which we obtain 

1 


C un = 




or 


arcsin(a) 

— e~ 2 ' 


—ia 1 
1 ia 


r) 




1 


sh(o:) 


y/2 s/d 1 - 1 \ ie 

1 


es 

— 77 arcosh(a) 


g— ^ arcosh(o;) 


VVct 2 - 1 \ ie2 

and for the metric operator 


arcosJ 


h(a) 


c br = 


sj a 2 — 1 


—1 a 
a 


( 21 ) 


( 22 ) 


Similarly, we have for a > 1 in the broken regime, 


(23) 


(24) 


Finally, we assume the system to be driven externally by 
varying the parameter a t . This is motivated by the optics 
experiment of Ref. [7], where one would change the optical 
pumping, i.e., vary the gain coefficient, 7 . Such a driving 
could be implemented in the set-up [7] by modulating the 
refraction index over the length of the waveguide. Projec¬ 
tive measurements could be realized by measureing ampli¬ 
tude and phase of the field at both ends of the waveguides. 

For broken as well as for unbroken PT-symmetry the 
gauge field A t ( 6 ) reads. 


A = 


ih ( 0 — i ri t \ 

2 ( a 2 t - 1 ) [ioct 0 ) 


(25) 


Equations (18)-(25) are all the ingredients necessary to 
compute the quantum work distribution ( 8 ) explicitly. 
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FIG. 2. (color online) Real part (red, solid line) and imaginary 
part (purple, dashed line) of the irreversible work (Wj rr ) (27) for 
the linear protocol (26) and ao = 3/2, that is a system starting in 
the regime of broken VT -symmetry. Parameters are f3 = 1, h = 
1, r = 1, and k — 1. Inset illustrates the violation of the Jarzynski 
equality (15) in real part (blue, solid line) and imaginary part 
(green, dashed line). 

For the sake of simplicity we further assume that a t is 
changed linearly from an initial value Q: () to a final value 
cti during time r, 

a t = a 0 + (a 1 -a 0 )t/T. (26) 

In Fig. 1 we plot the irreversible work, 

(WU = (W) - A F (27) 

as a function of the final value aq for a system that starts 
in the regime of unbroken VT- symmetry - namely a 0 = 
1/2. We observe that (Wi rr ) is always non-negative, as 
it should be as an expression of the second law of ther¬ 
modynamics, and that (W m ) diverges at the critical point, 
oq = 1. Similar behavior has been observed for quan¬ 
tum phase transitions [42], where the irreversible work 
can be interpreted as a quantum susceptibility. Diverging 
susceptibilities are a common feature of quantum critical 
points [43]. However, in quantum phase transitions the 
irreversible work usually exhibits a singularity, i.e., after 
passing through the transition (W- m ) drops back to a finite 
value [42, 43]. In the present case (W m ) diverges, sig¬ 
nifying a quantum critical point, but stays at infinity even 
after passing through the phase transition. In other words, 
the phase transition is thermodynamically “inhibited”, as 
the system’s response “freezes out”. The inset illustrates 
the validity of the quantum Jarzynski equality for systems 
with unbroken VT -symmetry (15) 

Figure 2 shows the same quantities, but for a system 
starting in the regime of broken VT -symmetry, here a 0 = 
3/2. The first observation is that (Wirr) has real and imag¬ 
inary parts. This is no surprise as the quantum work is de¬ 
fined as an average over differences of eigenenergies (7). In 


the broken regime, however, the eigenengeries are complex 
(20), and thus also (Wi rr ) has to be complex. This signi¬ 
fies the failure of the two-time measurement approach to 
describe the thermodynamics of systems with broken VT - 
symmetry. Nevertheless, we also observe that the abso¬ 
lute value of the irreversible work diverges at the critical 
point, oq = 1 , and that also the transition from broken into 
unbroken regime is thermodynamically inhibited. Finally, 
the inset illustrates the (complex) violation of the quantum 
Jarzynski equality, which is in full agreement with the gen¬ 
eral theory above - the quantum Jarsynki equality does not 
hold if the dynamics is not unital [34, 38-40]. 

Concluding remarks. In the present work we have 
shown how the quantum Jarzynski equality generalizes 
to VT -symmetric quantum mechanics. We have found 
that for quantum systems with unbroken symmetry the 
Jarzysnki equality holds, while this is not the case in the 
broken regime. The crucial requirement is that the dy¬ 
namics is at least unital. Hence, in time-dependent VT - 
symmetric quantum mechanics the Schrodinger equation 
has to be generalized including a gauge field so that all 
norms are preserved during the driving. 

These findings have been further analyzed for an ex¬ 
perimentally relevant system. We have illustrated that 
the VT -symmetric Jarzynski equality could be experimen¬ 
tally studied in optical set-ups consisting of two coupled 
waveguides. In these systems the time-dependent dynam¬ 
ics is frozen in the amplitude profile of the fields. Time- 
dependent driving could be implemented by modulation of 
the refraction index over the length of the waveguide. For 
this system we have found that the phase transition between 
the regimes of unbroken and broken VT -symmetry is ther¬ 
modynamically inhibited as the irreversible work diverges 
at the critical point - the system’s response “freezes out”. 

We emphasize that we chose the optical system due to 
its mathematical simplicity. However, all reported findings 
are completely general and also apply, for instance, to VT - 
symmetric systems in microwave billiards [44], photonic 
lattices [45], LRC circuits [46], optical lattices [47], meta¬ 
materials [48], or phonon lasers [49]. 
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